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The  h-p  Version  of  the  Boundary  Element  Method  with  Geometric  Mesh  on  Polygonal  Domains 

!.  Babuska1,  B.  Q.  Guo2  and  E.  P.  Stephan3 


Summary 

This  paper  applies  the  techniques  of  the  h-p  version  to  the  boundary  element  method  for 
boundary  value  problems  on  plane  non-smooth  domains  with  piecewisely  analytic  boundary  and  data. 
The  exponential  rate  cf  convergence  of  the  boundary  element  Galerkin  solution  is  obtained  when  a 
geometric  mesh  refinement  is  used  near  the  vertices. 

1.  Introduction 

The  h,  p  and  h-p  versions  (on  a  quasiuniform  mesh)  of  the  boundary  element  Galerkin  method 
for  integral  equations  on  polygons  has  been  studied  in  various  papers,  e.g.,  in  (1,  2.  3,  -4].  In  all 
versions  the  algebraic  rate  of  convergence  of  the  Galerkin  solution  is  restricted  by  the  vertex 
singularities  of  the  solution  of  the  integral  equations  although  it  could  be  very  smooth  away  from  the 
corners.  Based  upon  a  regularity  analysis  (ir>  countable  normed  spaces)  for  the  solution  of  the  integral 
equation  we  show  in  [5]  that  an  exponential  rate  of  convergence  with  respect  to  the  number  of  degrees 
of  freedom  can  be  achieved  for  the  h-p  version  by  simultaneously  reducing  the  mesh  size  and  increasing 
the  polynomial  decrees  of  the  boundary  elements  if  a  geometric  mesh  refinement  towards  the  vertices  is 
used.  Here  we  report  the  main  results  from  [5]. 

Let  C  be  a  bounded  polygonal  domain  with  vertices  A,,  1  <  i  <  M.  the  boundary  dP  be 
a  piecewise  analytic  curve 

M  _ 

op  =  r  =  u  r, 

i=l 

where  [V  is  an  open  line  segment  connecting  A,  and  A1  +  [  (A^^j  =  Aj).  By  u>,  we  denote  the  interior 
angle  at  A,. 

Let  k  >  0  integer,  denote  the  usual  Sobolev  spaces  furnished  with  the  norm 
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where  a  =  (a, .a-,),  a,  >  0  integer.  |q|  =  a,  +  a2  and  DQu  =  — °Q  UQ?.  The  space 
defined  as  the  restriction  of  functions  in  H*(Q),  i.e.. 


,1-1/2, 


dx1  <9x 


Hfc-1/2(r>  =  {u|r.  u  eHfc(ft)}.  for  k  > 


H*  1/0( r )  =  L2(F),  for  k  =  i. 


H*  Xj  "(T)  =  (H  ‘ 7 “  f ) ) ^  (dual  spaces),  for  k  <  5. 

For  the  investigation  of  singularities  at  corners  we  introduce  weighted  Sobolev  spaces  and  countable 

normed  spaces  on  the  hmmd-’y  r  r  "*  !  —  (a,b),  and  fo;  A  c  (a.b),  ij  _  jA  —  aj,  i2  .=  ]A  —  b|. 

,  3 .  + 1  -  -  -  - 

<J>-  Ax)  =  Jj  r  .(x),  3  =  (3i,3->).  0  <  3u3->  <  1 ,  k  integer.  Now  we  define  for  k  >  £  >  0.  and 
i  =  l 

integer  £  >  0: 

H*  £d)  =  {u  6  H£-‘(I)  if  £  >  0:  !!<i>3  +  m_£(x)u(m,(x)||L2(1)  <  00.  for  £  <  m  <  k} 

B£(l)  =  {u  €  V  k  >  £,  |!4.+t_£u(t,(x)||  <  Cdl“£(k-e).  C  >  1,  d  >  1  independent  of  k}. 

k.e  i  .  .  .  £ 

For  any  T,  C  F.  Hj  (T,)  and  B^  (TJ  are  defined  correspondingly  with  3,  =  {3,\A3l2)-  B^(T)  = 

Mr  . 

n  B . * (r, )  With  £  =  (£,.£2 . £m).  3  -  (3i . ^yj).  We  shall  write  3,  >  3,  (resp.  3  >  3)  if 

1=1  .  i,U\ 

3,  }  >  3,  },  j  =  1,2  (resp.  3,  >  3,,  1  <  i  <  M).  By  B^  (T)  we  denote  for  1  <  i  <  M  the  space 
£  £  +  1 

n  B3,(r.)  x  n  b*.  (r,). 

0<3,<1  i<3,<l 

t,£  v  i  >  — 

Remark  1.1.  Hi  (T)  an  !  1  (T)  are  the  trace  spaces  of  functions  belonging  to  the  weighted  Sobolev 

k,t  k,i 

space  (Q)  and  the  countable  normed  space  (Q),  respectively  (see  [6.  7,  8]) 


In  the  next  sections  we  analyze  the  regularity  of  the  solution  of  integral  equations  for  Neumann 

£,£+ 1  ;a- 

and  Dirichlet  boundary  value  problems  in  terms  of  countable  normed  spaces  B^  (T),  design  the 
geometric  mesh  and  the  distribution  of  the  degrees  of  polynomials  in  the  boundary  element  Galerkin 
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method,  which  lead  to  the  exponential  rate  of  convergence  with  respect  to  the  number  of  degree  of 
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freedom.  For  the  proof  of  the  theorems  we  refer  to  [5]. 


The  Neumann  Boundary  Value  Problem 


We  ccnsiJei  the  Neumann  problem  on  a  polygonal  domain  Q 


I  u  =  0  in  Q 
<9u| 

T 


(2.1) 


dn  i  ® 


where  ^  means  the  normal  derivative  with  respect  to  the  outer  unit  normal,  and  e  satisfies 
cm 


J  gds  =  0 

r 

Let  D  and  be  the  integral  operators  for  x  £  T 

Dufx)  =  -  f  £;  l  ^  («n|x  — y|)u(y)ds(y\ 
and  F 


K'$j|(x)  =  ?&_  f  ^Kx-yi 


<9u(y) 

0nv 


My) 


VHZ  p 

Then  the  first  kind  integral  equation  for  the  Neumann  boundary  value  problem  (2.1)  reads  as 


(2.2) 


Du  =  f  on  T 


(2-3) 


with  f  =  (1  —  K?)g,  see  [1]. 

0,1 

Theorem  2.1.  [5]  For  given  g  €  (T)  satisfying  (2.2),  the  integral  equation  (2.3)  together  with  the 

12 

side  condition  J  uds  =  0  has  a  unique  solution  u  E  Bj  (T),  for  some  0,  0  <  0  <  1  where  0  depends 

r 

on  3  as  well  as  on  the  geometry. 

Now  we  discuss  the  numerical  solution  of  (2.3)  by  the  h-p  version  of  the  boundary  element 
Galerkin  method  with  a  geometric  mesh. 

Let  ft  be  a  L-shape  domain  shown  in  Fig.  2.1.  We  assume  for  simplicity  that  the  solution  u  of 
(2.3)  belongs  to  b|(T)  (resp.  B^(T))  with  =  |y|Pl’2,  4>-^  =  |x|*2,1,  |  j  <  1  (resP-  g  < 

0 1 2<02  i  <  a)  and  $  a  =  1  for  3  <  j  <  6,  i.  ,  the  singularity  occurs  only  at  the  origin.  For  example 
i.  P  i 
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2/3  0  i/3  fl 

this  is  the  case  of  u  =  t  cos  =  6  (resp.  u  =  r  sin  on  r  where  (r,8)  denote  the  polar 
coordinates  centered  at  the  origin  (see  [9,  lOj). 

Let  cr  £  (0.1)  be  the  mesh  factor  and  n,  integer,  be  the  number  of  layers,  and  let  T,  1  <  i  < 
I  ( j  J .  i  <;  j  <  n+1  be  the  boundary  internals  such  that  dist(0,ri;)  -  crn  +  1_;,  1  <  j  <  n  +  1  and 
dist(0.r,  j)  =  0,  1  <  i  <  I(j).  Then  T  a  =  {T,  },  1  <  i  <  I  (j ) ,  1  <  j  <  n  +  1}  is  called  the  geometric 
mesh  on  associated  with  <7  and  n.  Fig.  2.2  shows  a  sequence  of  the  geoemtric  meshes  with  cr  ~  0.15. 

Let  P  =  {p,  1  <  i  <  I(J1.  1  <  j  <  n+1}  be  the  degree  vector  with  p,  >  1  integer.  The 

boundary  element  space  associated  with  the  geometric  mesh  r„  and  degree  vector  P  is  defined  by 

P  n 

S  (T-l  =  {<p|<p|p  is  a  polynomial  of  degree  <  p,  } 

1  i.;  ' 

and 

sp(r")  =  sp(r")  n  c°(r)  c  B1/2(n 

The  boundary  element  Galerkin  procedure  for  the  integral  equation  (2.3)  with  given 

0  I  1  O  p  fl 

g  €  Bj  (T)  is  to  find  up  €  S  (T,)  satisfying  J  upds  =  0  such  that 

r 

(Dup.w)  =  ((1-K')g,w),  V  w  6  SP(r”)  (2.4) 

where  {•,■)  denotes  the  duality  between  H  1/2(T)  and  li1/2(r).  \V  e  have  the  following  approximation 
result  of  the  h-p  version  of  the  boundary  element  Galerkin  method. 

2  1 

Theorem  2.2.  [5]  Let  u  G  B^(I')  (resp.  Bg(T))  be  the  solution  of  the  integral  equation  (2.3)  and  T  be 
the  boundary  of  the  L-shaped  domain  shown  in  Fig.  2.1,  0,  -  =  01  j  =  /32  2  =  0,  3  <  i  <  6,  j  =  1.2. 
g  <  d]  2,d21  <  1,  (resp.  g  <  P\  2'  ~02  \  <  5)  and  Ter.  (7  €  (0,1)  be  the  geometric  mesh  on  T.  Lei 
SA  (T,J  denote  the  boundary  element  space  defined  above  with  pi;  =  p;  >  1,  j/i  <  p;  <  v n,  0  <  p  < 
u  <  oo.  Then  the  boundary  element  Galerkin  solution  <  6  SP(r")  of  (2.4)  converges  to  u  in  HI/2(T) 

exponentially,  i.e., 

||u-up||  1/2  <  Ce“bNl/2 

H  7  (T)  _ 


where  N  is  the  number  of  degrees  of  freedom,  C  and  b  are  some  constants  independent  of  N. 
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3.  The  Dirichlet  Boundary  Value  Problem 


In  this  section  we  consider  the  Dirichlet  boundary  value  problem 


Au  =  0  in  fi 
u|r>  -  2 


With  the  integral  operators 


V  ~  f  I  75"  «n|x-y|ds(y),  x  €  fi 


P  dny 


Ku(x)  =  -  |  J  A  {fnjx— y|}u(y)ds(y),  x  £  Q 


(3.1)  leads,  as  shown  in  [2],  to  the  first  kind  of  integral  equation 


v  &  =  f 


with  f  =  (1  +  K)g,  for  which  there  holds  the  following  result 

Theorem  3.1.  [5]  Let  cap(f)  ^  1  where  cap(T)  is  the  analytic  capacity  of  T.  Then  for  given 

g  6  Bj  (F)  O  C°(F)  there  exists  exactly  one  solution  ^  £  B.j  (F)  of  { 3.2)  with  some  3,  0  <  3  <  1. 

Now  we  consider  the  Galerkin  solution  of  (3.2)  obtained  by  the  h-p  version.  For  simplicity  we 
assume  again  that  Q  is  the  L-shaped  domain  shown  in  Fig.  2.1  and  the  solution  of  (3.2)  has  a 

singularity  at  the  origin  only.  Then  the  geometric  mesh  on  T  and  the  boundary  element  space 

SP( )  are  defined  as  in  the  previous  section.  Obviously,  S^>(r£T)  C  L2(T)  C  H  1/2(T). 

P  n 

The  Galerkin  procedure  for  the  integral  equation  (3.2)  is  to  seek  ipp  £  S  (r<,)  such  that  for  all 

wP  £  sp(r") 

(ViPp,<t>p)  =  ((1  +  K)g,*p)  (3.3) 

For  boundary  element  solutions  ifp  we  have  the  following  approximation  theorem. 

Theorem  3.2.  [5]  Lei  ^  £  B^(T)  ( resp .  B^(T))  be  the  solution  of  the  integral  equation  (3.2)  and  T  is 
the  boundary  of  the  L-shaped  domain  shown  in  Fig.  2.1,  with  cap(T)  ^  1,  and  01,1  —  02,2  ~  0i,i  for 

3  <  «  <  6,  ;  =  1,2,  jj  <  0 1,2'0 2,i  <  1.  ( resp .  |  <  0i,2i02,i  <  \)  and  lei  r*>  a  e  C0’1)  ie  tflc 
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P  n 

geometric  mesh  on  T  and  S  (r„)  be  the  boundary  element  space  defined  above  mth  p,  —  P;  >  0. 

P  O 

JP  £  P;  £  "t1  0  <  p  <  t'  <  x;.  Then  the  boundary  element  Galerkin  solution  vp  (E  S  (T,)  of  (3.3) 
converges  to  ^  in  H  '  ‘(f)  exponentially,  i.e.. 


H,,-  - 


<  Ce~bN‘/’ 
on  H  /'(T)  “ 


uhere  N  is  the  number  of  degrees  of  freedom.  C  and  b  are  some  constants  independent  o/S. 


•4.  Conclusion 

The  regularity  results  for  the  solutions  of  the  integral  equations  for  the  Dirichlet  and  Neumann 
boundary  value  problems  of  the  Laplacian  can  be  generalized  to  bipotential  and  elasticity  problems 
with  essential,  natural,  and  mixed  boundary  conditions.  The  h-p  version  of  the  boundary  element 
method  possesses  advantages  over  the  finite  element  method  such  as  reducing  the  number  of  degree  of 
freedom  and  avoiding  the  difficulties  in  the  treatment  of  non-homogeneous  essential  boundary- 
conditions  in  the  finite  element  method  (see  [10]).  Although  the  geometric  mesh  show-'  in  Fig.  2.1  is 
designed  for  the  problems  with  singularity  at  one  corner  it  can  be  generalized  without  any  difficulty  to 
the  case  that  the  singularity  occurs  at  several  corners  of  T,  and  the  exponential  rate  of  convergence  can 
be  proven  as  well.  All  theorems  above  will  hold  if  Q  is  a  curvilinear  polygon. 
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